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Abstract 

We obtain an exact result for the soft gluon production and its pt distribution due to a gluon 
loop in a constant chromo-electric background field E a with arbitrary color. Unlike Schwinger's 
result for e + e~ pair production in QED which depends only on one gauge invariant quantity, the 
electric field E, we find that the px distribution of the gluons depend on two gauge invariant 
quantities, E a E a and [d abc E a E b E c f . 
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Although hard gluon production at high energy colliders is computed by using perturba- 
tive quantum chromodynamics (pQCD), it is not possible to calculate soft gluon production 
by using perturbation theory because the coupling constant becomes large at low energy. 
In this paper we study non-perturbative soft gluon production due to vacuum polarization 
in a constant chromo-electric field by applying the background field method to QCD with 
gauge group SU(3). While e + e~ production from a constant electromagnetic field in spinor 
QED has been calculated long ago by Euler and Heisenberg [l], and Schwinger 2!\, and by 
Weisskopf for scalar QED, there is of course no direct photon production in QED as 
photons do not interact with a classical electromagnetic field. In the case of QCD, gluons do 
interact with a background chromofield and gluons are produced. As we want to study the 
Pt distribution of gluons produced at collider experiments, we shall not employ Schwinger's 
proper time method but, rather, we will directly evaluate the path integral and obtain the 
partition function which gives both the px distribution and the total gluon production. 

For the interpretation of experiments which probe the quark-gluon plasma at high-energy 
large-hadron colliders such as RHIC (Au-Au collisions at y/s = 200 GeV) [J and LHC (Pb-Pb 
collisions at ^/s = 5.5 TeV) [5J] it might be necessary to know the pt spectrum of soft gluons 
produced by a chromo field. The physical picture behind this calculation is that two heavy 
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nuclei collide and then move apart, creating a classical chromo field in between 
consider here a constant chromo-electric field along the beam direction. Subsequently this 
constant chromo-electric field breaks up into quark_anti-quark pairs and gluons. The quark 
anti-quark production has been calculated be 
has been approximated by WKB methods 



'ore p and the corresponding p T distribution 
Also the total gluon production rate by a 



covariantly constant field (dN/d x) has been calculated for SU(2) in [10j,|ll| and for SU(3) 
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12j . The pr distribution (dN/d xd pr) of gluon production due to a constant chromo- 
electric field which we present in this article is new, as far as we know, and we consider the 
realistic case of SU(3). 

We obtain the following formula for the number of non-perturbative soft gluons produced 
per unit time and per unit volume and per unit transverse momentum from a given constant 



chromo-electric field E a 



This result is gauge invariant because |Ai|, |A 2 |, and |A 3 | are the positive square root of the 



following gauge invariant real positive quantities 
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A 2 = y [l + cos V3 

= y[l + cos(| + e)], (2) 

where 9 is real and is given by 

cos 3 e = -l + 6C 2 /cf . (3) 

They depend only on the Casimir invariants for SU(3) 

Ci = E a E a , C 2 = [d abc E a E b Ef, (4) 

where a, b, c = 1,...,8 are the color indices of the adjoint representation of the gauge group 
SU(3). Note that 9 is real because Cf - 3C 2 > 0. 

This can be contrasted with the corresponding formula for massless fermion pair produc- 
tion 

dN e+e - -\eE\ , , 



ln[l - e"^]. (5) 



dtd 3 xd 2 pT 47r 3 

The result for the pr distribution in @ was obtained by WKB methods in|9| and integration 
over pt reproduces Schwinger's result for total production rate dN/ d 4 x j2j . (The first term 
in the expansion of the logarithm of the total production rate dN/ d 4 x was already obtained 
by Heisenberg and Euler yj. For an excellent review on Euler-Heisenberg actions, see 
In our result the symmetric tensor d a b c appears. Hence the extension of Schwinger's formula 
for e + e~ pair production to gluon production is not straight forward. We now present a 
derivation. 

The Lagrangian density in the background field method of QCD with covariant back- 
ground gauge fixing term in Feynman-'t Hooft gauge is given by |14] 

C gluon = \ Q» a M*[A] Q» b (6) 



where 



M£[A] = ^ u [D p (A)D"(A)} ab - 2gt bc F; v (7) 
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with rj^ = (—1, +1, +1, +1) the Minkowski metric. The corresponding ghost Lagrangian 
density is given by 



C gh ost = TDf[A]D^[A + Q] x c = TK ab [A,Q] x h 



(8) 



where Q, A and x are the gluon quantum field, the background field, and the ghost field, 
respectively, and D^[A]Q« = [d a Q a u + g f abc A\Ql\. 

The one-loop effective action for a gluon loop in a background field A® is given by 



$L = - i\n(Det M)~ 1/2 = % - Tr[\nM[A\ H - lnM[0] H] 



(9) 



where H^ u = rf . We added the matrices H^ v since this will allow us to factorize the trace 
over Lorentz indices. The trace Tr contains an integration over d A x and a sum over color 
and Lorentz indices. We replace the logarithm by using the exponential representation 

ds , 



= - Tr 

gluon 2 



f°° as 
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e is (M[0] H+ie) _ e is (M[A] H+ieh 



(10) 



We shall say more about the s-integration contour later. We assume that the electric field is 
along the z-axis (the beam direction) and we choose the gauge Aq = so that A3 = —E a x°. 
The color indices (a=l,....8) are arbitrary. Since A ab = if abc E c is hermitian and anti- 
symmetric, its eigenvalues are real and come in pairs (A, — A). We will show later that the 
matrix A ab has six non vanishing eigenvalues. So after diagonalization it reads 



(11) 



We replace the derivative d n by ip u satisfying [p u ,x u ] = —iS 1 ^ and obtain 



M[A\*rT = S; [(po 2 - PT 2 ~ Ps 2 )S ab - 2gA ab x°p 3 - g 2 (A 2 ) ab x° 2 ] + 2igA ab S " (12) 



where pt = \jp\ + p\ is the transverse momentum of the gluons (transverse to the electric 
field direction) and 
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(13) 
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To reduce this problem to one harmonic oscillator, we make the usual similarity transforma- 
tion [151] . We also make a similarity transformation in group space which diagonalizes the 
matrices A 

M[A]%rT = {e ip3po ^ [(pV-PT 2 ~ g\^)x° 2 ) 5; + 2igA d S;\ e'^^f . (14) 

Recalling that A ab has two vanishing eigenvalues, we should read this expression as follows: 
for a=b=l,..6 each denotes an eigenvalue ±Ai, ± A2 ± A3, while for a=b=7,8 one finds 
an expression which is independent of A a ^ and which we omit because it cancels against the 
contributions from M[0]. 

The trace over Lorentz indices factorizes because 5" commutes with S,, u and we obtain 



[2 + 2cosh2sgX j ] - Ae is{p '° 2 - p ^ +ie) ] (15) 

The trace tr denotes integral over a complete set of x eigenstates. From here we proceed as 
in QED [15]. We add complete sets of p eigenstates, and obtain 

00 ds 1 f J4 _ f J2 „ _-i«2, — set 1 _ > , 1 + cosh2sg\ J 



tlS ■ r ip 3 Po/g^j p is(P0 2 -PT 2 - g 2 \ 2 ,x 2 +ie) -ip 3 po/gX 3 



The one-loop effective action for the ghost in the background field is given by 

SSL* = -MDetK) = -iTr —\e is - e is (17) 

y Jo s 



where i^ ab [A] is given by (JHJ). Since there is now no trace over /1, v but also no square root of 



abc £pc 
flu 



the ghost determinant, we find an overall factor —1/2. In addition there is no term f abc F^ 
in K[A] so that the term cosh2sg\j is absent. We obtain then 

41< = 5 f t £ (SF I d ' x l*» ~ "^sssj^i - (18) 

Adding the effective action due to a gluon and a ghost loop we find for the one- loop effective 
action 

2 Jo s ^ (2vr) 3 7 J uy Jl sin/w»|</Aj| s J V ; 

The s- integral at fixed px is convergent at s — > 0, but integration over pt yields an extra 
factor 1/s, and charge renormalization cures this ultraviolet problem by substracting also 



the term 



function 10, [ll ] . 



inear in s in the expansion of cosfc2 ^ A J I _ j n f ac ^ the term linear in s yields the beta 

c sinhs\g\j \ J 



W gluon = 2ImSU = ' J d 2 p T ]T (-l)^ 1 ± \£M e [ -W ] . (21) 



(27T) 



A more serious problem occurs for s — > oo. Obviously, the integral over s diverges, 
indicating the presence of infrared divergences. Although such divergences in Feynman 
graphs with on-shell external gluons are well known, the full non-perturbative treatment of 
these divergences seems complicated. However, we are only interested in the imaginary part 
of the effective action, since it gives us the gluon production rate, and hence we proceed 
similarly to Yildiz and Cox (see eq. 22 in [l2|). We use the well-known expansion 

1 00 (—l) n 

1/sinhx = - + 2x Y V^7 (20) 

x ~l x l + n z n z 

and then we formally replace s by —is (as first advocated by Schwinger in QED where no 

infrared problems are present). The integral is now real, except for half-circles around the 

poles at slgAjl = — irm for n=l,2,3... The term 1/x in ()20|) cancels against the term 1/s 

in (fTTJj) . This yields probability for gluon production per unit time and per unit volume 

n=l j=l 

All that is left is to determine the eigenvalues Xj (j = 1,...,8) of the matrix if abc E c . First 
we evaluate the determinant of the 8x8 matrix \J ahc E c — A 5 ab ] and find 

Det[f abc E c - A 5 ab ] = A 2 [A 6 + AX 4 + BX 2 + C] 

= A 2 (A 2 + A 2 ) (A 2 + A 2 ) (A 2 + A 2 ). (22) 

Because E b is an eigenvector of the matrix if abc E c , and because eigenvalues come in pairs 
(Xj, —Xj), there are two eigenvalues zero, which explains the factor A 2 . The gauge invariant 
quantities A, B, and C in the above equation can only depend on the Casimir invariants in 
(JJJ) and we find 

A = ^E a E a , B = 4-, C = 4 ([E a E a ] 3 - 3 [d abc E a E b Ef), (23) 
z 4 lb 

where d a b c is the symmetric invariant tensor in the adjoint representations of SU(3). It 

follows that the three eigenvalues satisfy 

A? + A 2 + A 2 = A 

AiA 2 , + A2A 2 + A 2 A 2 = B 

A 2 A 2 A 2 = C, (24) 



the solution of which is given by eq. (0). 

In this letter we have derived the pr distribution of the soft gluons produced by the 
vacuum polarization due to a gluon loop in a constant chromo-electric field. We have used 
the background field method of QCD with gauge group SU(3). For application at RHIC 
and LHC we have constructed the pt distribution of the gluon production. We find that, 
unlike the case in QED where the e + e~ pair production rate depends on one gauge invariant 
quantity the pr distribution of the gluon production rate depends on two gauge invariant 
quantities, E a E a and [d a b c E a E b E c ] 2 . We intend to add a chromo- magnetic field later. 
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